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Digital Computer Laboratory- 
Massachusetts Institute of Technology 
Cambridge, Massachusetts 

SUBJECT* GROUP 61 SUBROUTINE LIBRARY 

Tot Group 61 

From* W. Lone 

Date* October 2, 1952 

Abstract* Additional subroutines are given and should be appended to M-16S1, 
dated September 5, 1952. Corrections and clarifications of that 
memorandum are presented. 

The following subroutines are included* 

SRL #9 4x7 Decimal Display 

#10 Radius Vector & Azimuth Subroutines 
#11 Unscrambling a Word from Magnetic Tape 
#12 Arcsine in Revolutions 
. #13 Sine-Cosine Subroutine 
SRL #4 ARCSINE COMPUTATION has been corrected as follows* 

106 ad 137 
137 0.200000 

It is pointed out that since 9, which lies between - -k and ♦ «•- , has a scale 
factor of Z" f arcsin x lies between - •»• and + -s- • 

In SRL #3 the constant 0.11404 in register 202 is ■% • 

The following corrected subroutine is included* SRL #4 Arcsine 
Computation. 



Signed* //. <£<n^£ 



W. Lone 



frttiJu 



Approved * k ' 'l v C/y U/Qjj\ 
C. R. Wieser 



WL/efh 

Attached forms* 

SRL #9 4x7 Decimal Display 

#10 Radius Vector & Azimuth Subroutines 
#11 Unscrambling a Word from Magnetic Tape 
#12 Arcsine in Revolutions 
#4 Arcsine Computation 
#13 Sine-Cosine Subroutine 



Title » SHL #9 4 x 7 DECIMAL DISPLAY 

Abstract t Enter with In| x 10" in AC, where k a number of digits to be 
displayed. Ho sign is displayed. 



100 
101 
102 
103 
104 
105 
106 
107 
110 
111 
112 
113 
114 
115 
116 
117 
120 
121 
122 
123 
124 
125 
126 
127 
130 
131 
132 
133 
134 
135 
136 
137 
140 
141 
142 
143 
144 
145 
146 
147 
150 
151 
152 
153 
154 
155 



ta 150 

ad 232 

ts It 

ca 224 

ts 3t 

ca It 

mh 230 

ad 233 

td 115 

si 17 

ts It 

sp 114 

ta 144 
ca— — 

sp 117 

ta 140 

°GS MU 

cs 227 

ts 4t 

ca 3t 

si 11 

si 600 

ca 2t 

sf 5t 
su 231 

ts 2t 

ca 5t 

ad 4t 

cp 151 

ao 3t 

ca 2t 

si 1 

sp 

ca 115 

ad 230 

td 115 

sp 

ao 3t 
su 225 

cp 105 

sp 

ts 4t 

ad 226 

si 11 

re 

sp 126 



200 


0.37301 


201 


0.00302 


202 


0.74711 


203 


0.40711 


204 


0.07610 


205 


0.43705 


206 


0.37311 


207 


0.00601 


210 


0.33311 


211 


0.23311 


212 


0.40676 


213 


0.77700 


214 


0.44706 


215 


0.44666 


216 


0.04377 


217 


0.42671 


220 


0.44662 


221 


0.76603 


222 


0.44666 


223 


0.44676 


224 


P x 

r 


225: 


P * f 


226 


py o 


227 


p 6 


230 


p 10 


231 


.40000 


232 


1/2 x 10 


233 


si 200 



abscissa of desired left 
extremity of display 
abscissa of desired right 
extremity of display 
ordinate of upper 
extremity of display 



-k 



remaining to be 



It digits of 
displayed 
2t remainder of constant 
3t present x deflection 
4t y counter 
5t sf 



Titles 



SRL #10 RADIUS VECTOR & AZIMUTH SUBROUTIHES 



Abstract* Enter with «■ in 200, 2. i n 201. The output will be x in 2t, 

(p in It where -7p ~ (b 6l -^ revs . Should radius vector alone be 

desired the changes are 100 ta 120; 120 sp with orders 
beyond 120 not needed* 



Formulas 



E = /x 2 *y 2 = »(|x| ♦ |y| - b jj^'fjrl ) 



< lj = sgn x 



fr .. *■■«"> hfffe .] 



Precisions in R, 0.75$; in(b , max. error 1.23 degree. 



100 


ta 


136 


101 


cm 


200 


102 


ta 


It 


103 


cm 


201 


104 


ad 


It 


105 


ts 


2t 


106 


cm 


201 


107 


mr 


It 


110 


mr 


202 


■111 


dv 


2t 


112 


si 


17 


113 


ts 


It 


114 


si 


1 


115 


ad 


2t 


116 


mr 


203 


117 


ts 


2t 


120 


cs 


201 


121 


cp 


124 


122 


ad 


It 


123 


sp 


125 


124 


su 


It 


125 


mr 


204 


126 


dv 


2t 


127 


si 


15 


130 


ad 


205 


131 


ts 


It 


132 


es 


200 


133 


ep 


136 


134 


cs 


It 


135 


ts 


It 


136 


sp- 





200 ~ 



201 



2 



202 1.32404 -b/2 

203 0.40172 a/2 

204 0.37176 c/2 

205 0.17602 d/2 



Title: SRL #11 UNSCRAMBLING A TflORD FROM MAGNETIC TAPS 

Abstracts The word to be unscrambled is in AC, and at the conclusion the 
unscrambled word is in AC. BR must be clear when entering this 
subroutine. This subroutine takes about 2.1 milliseconds. The 
maximum permissible time is 2.5 milliseconds. 



100 


ta 145 


101 


clh 22 


102 


ts t- 


103 


su t- 


104 


clh 2 


105 


tst 2 


106 


sut 2 


107 


clh 2 


110 i 


ts t 3 


111 


SU V 


112 


clh 2 


113 


ts t 4 


114 


su t 4 


115 


clh 2 


116 


ts t g 


117 


su t c 
5 


120 ; 


clh 2 


121 


"CS t g. 


122 


SU tg 


123 


clh 2 


124 


ts t„ 


125 


SU t„ 


126 


clc 2 


127 


sic 2 



130 


ad t- 


131 


sic 2 


132 


adt 6 


133 


sic 2 


134 


ad t_ 
5 


135 


sic 2 


136 


ad t 4 


137 


sic 2 


140 


ad t„ 


141 


sic 2 


142 


ad tg 


143 


clc 2 


144 


ad t.. 


145 


sp 



Title: SEL #12 ARCSINE IN REVOLUTIONS 

Abstract* Enter -with x in AC (-1<£ x < l). The output of 9 = arcsin x 
in revolutions is in AC where (-1/4^ ^< l/4). The maximum 
error is about .001 rev^.5 degree* 



100 ta 131 

101 ts It 

102 cm It 

103 ts 2t 

104 cs 2t 

105 sr 2 

106 ad 204 

107 sp 100 (SRL #7) 

110 ts 3t 

111 ca 2t 

112 mh 203 

113 sr 1 

114 ad 202 

115 mh 2t 

116 sr 2 

117 ad 201 

120 mh 2t 

121 sr 3 

122 ad 200 

123 mr 3t 

124 su 204 

125 ts 3t 

126 ca It 

127 cp 132 

130 cm 3t - 

131 sp 

132 ca 3t 

133 sp 131 



200 


0.40000 


201 


1.35330 


202 


0.55340 


203 


1.50100 


204 


0.20000 



Titles SRL #4 ARCSINB COMPUTATION 



Abstracts Enter with x in AC, Q& x < 1. The output is 9 » 2 (sin x) 

radians in AC (— g- Sdi?). It follows then that - -k £ sin" x ^ 
+ jr • This must be used in conjunction with "Square Root Sub- 
program > w SEL #7. 





Precisions 


100 


ta 127 


101 


ts. It 


102 


cm It 


103 


ts ; 2t 


104 


cs >2t 


105 


sr 2 


106 


ad 137 


107 


sp 100 (SEL #7) 


110 


ts 3t 


111 


ca 2t 


112 


mr 135 


113 


ad 134 


114 


mr 2t 


115 


ad 133 


116 


mr 2t 


117 


sr 1 


120 


ad 132 


121 


mr 3t 


122 


su 136 ... 


123 


ts 3t "' 


124 


ca It 


125 


cp 130 



Maximum error about #001 radian* 



126 
127 
130 
131 
132 
133 
134 
135 
136 
137 



cm 3t 

sp 

ca St 

sp 127 

+.7853 

-•2121 

+.0743 

-.0187 

+.3927 

0.20000 



Title $ 



SEL #13 SUSE-COSINE SUBROUTINE 



Abstracts This variation of SEL #1 enables one to enter the subroutine with 
© such that - l/2 £ < + l/2 or £ © < +1 measured from Bast. 
Q should be in AC when entering the routine. Output will be 
sin © in 200, cos © in 201. 



100 


ta 125 


101 


sp 102 


102 


ta 123 


103 


cl 1 


104 


ts 200 


105 


cm 200 


106 


ts t x 


107 


mr 202 


110 


■ su 203 


111 


mr t-i 
ad 204 


112 


113 


mr ti 
ad 205 


114 


115 


mh 200 


116 


si 3 


117 


ex 201 


120 


ex 200 


121 


cl 31 


122 


sa 206 


123 


sp 


124 


si 


125 


sp 



200 


sin © 


201 


sin ©, cos 


202 


0.34540 


203 


0.71301 


204 


0.03671 


205 


0.30647 


206 


0.20000 



(clears spec, add mem.) 



